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Discrete Painleve equations and discrete KdV 
equation over finite fields 

Masataka Kanki* and Jun Mada^ and Tetsuji TOKIHIRO^ 

April 19, 2013 

Abstract 

We investigate some of the discrete Painleve equations (dPu , gFj and 
gPjj ) and the discrete KdV equation over finite field^. The first part 
concerns the discrete Painleve equations. We review some of the ideas 
introduced in our previous papers [U [2] and give some detailed discus- 
sions. We first show that they are well defined by extending the domain 
according to the theory of the space of initial conditions. We then extend 
them to the field of p-adic numbers and observe that they have a property 
that is called an 'almost good reduction' of dynamical systems over finite 
fields. We can use this property, which can be interpreted as an arith- 
metic analogue of singularity confinement, to avoid the indeterminacy of 
the equations over finite fields and to obtain special solutions from those 
defined originally over fields of characteristic zero. In the second part we 
study the discrete KdV equation. We review the discussions in [3] and 
present a way to resolve the indeterminacy of the equation by treating it 
over a field of rational functions instead of the finite field itself. Explicit 
forms of soliton solutions and their periods over finite fields are obtained. 

1 Introduction 

In this article, we study several discrete integrable equations over finite fields. 
One of the problems we encounter when we treat a discrete equation over a 
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finite field is tliat its time evolution is not always well defined. This problem 
cannot be solved even if we extend the domain from F^™ (p is a prime number) 
to the projective space PF^™ :— F^™ U {oo}, because PFpm. is no longer a field; 
we cannot determine the values such as ^^ O-oo, 00 + c» and so on. In this paper 
we study two strategies to define the time evolution over a finite field without 
inconsistencies: [I] One is to reduce the domain so that the time evolution will 
not pass the indeterminate states. [II] The other is to extend the domain so 
that it can include all the orbits. We take the latter strategy [II] and adopt 
two approaches: (i) The first approach is the application of the theory of space 
of initial conditions developed by Okamoto [U [5l El [7] and Sakai [8] . (ii) The 
second approach is to extend the domain of initial conditions to the field of 
j^adic numbers Qp. 

The first part of this article concerns the one-dimensional dynamical systems, 
in particular the discrete Painleve equations. Through the approach (i), we 
show that the dynamics of those equations can be well defined in the space 
of initial conditions even over the finite fields [5]. The approach (ii) is closely 
related to the so called arithmetic dynamics, which concerns the dynamics over 
arithmetic sets such as Z or Q or a number field that is of number theoretic 
interest [9]. In arithmetic dynamics, the change of dynamical properties of 
polynomial or rational mappings give significant information when reducing 
them modulo prime numbers. The mapping is said to have good reduction if, 
roughly speaking, the reduction commutes with the mapping itself [5]. The 
action of the projective linear group PGL2 give typical examples of mappings 
with a good reduction. The QRT mappings [10] and several bi-rational mappings 
over finite fields have been investigated in terms of integrability by choosing the 
parameter values so that indeterminate points are avoided [HI [121 [13] . They also 
have a good reduction over finite fields. We prove that, although the integrable 
mappings generally do not have a good reduction modulo a prime, they do 
have an almost good reduction, which is a generalised notion of good reduction. 
We first treat in detail the discrete Painleve II equation (dPn ) [H] over finite 
fields, on which we have briefly reported in our previous letter [1], and then 
apply the method to the g-discrete Painleve equations (gPi and gPn ). The time 
evolution of the discrete Painleve equations can be well defined generically, via 
the reduction from a local field Qp to a finite field ¥p. The theory is then used 
to obtain some special solutions directly from those over fields of characteristic 
zero such as Q or M. 

The second part concerns the two-dimensional evolution equation, in partic- 
ular, the discrete KdV equation (dKdV). Since the dKdV equation evolves as a 
two-dimensional lattice, the number of singular patterns we obtain is too large 
to be investigated even for a system with small size. Therefore, it is difficult to 
determine the space of initial conditions through the approach (i) above. In- 
stead, we take the second approach (ii). We can well-define the dKdV equation 
by considering it over the local field Qp just like the discrete Painleve cases. We 
introduce another extension of the space of initial conditions. Since the dKdV 
equation has a parameter 5, we can see S as an indeterminant (variable) and 
can define the dKdV equation over the field of rational functions F^ {6) where 



r = p™ (to > 1) [S'. We define the soliton solutions over tlie finite fields and 
discuss their periodicity. 

2 The dPii equation and its space of initial con- 
ditions 

The discrete Painleve equations are non-autonomous, integrable mappings which 
tend to some continuous Painleve equations for appropriate choices of the con- 
tinuous limit [13]. They are non-autonomous and nonlinear second order ordi- 
nary difference equations with several parameters. When they are defined over 
a finite field, the dependent variable takes only a finite number of values and 
their time evolution will attain an indeterminate state in many cases for generic 
values of the parameters and initial conditions. For example, the dPn equation 
is defined as 

M„+l + Uri-1 = " " „ {n £ Z), (1) 

where Zn ~ 6n + zq and a,(5, zq are constant parameters [T5]. Let r — p™ for 
a prime p and a positive integer to, G Z+. When ^ is defined over a finite 
field ¥r, the dependent variable u„ will eventually take values ±1 for generic 
parameters and initial values (uo,ui) G F^, and we cannot proceed to evolve it. 
If we extend the domain from F^ to (PFr)^ = (F^ U {oo})^, PF^ is not a field and 
we cannot define arithmetic operation in ([T]). To determine its time evolution 
consistently, we have two choices: One is to restrict the parameters and the 
initial values to a smaller domain so that the singularities do not appear. The 
other is to extend the domain on which the equation is defined. In this article, 
we will adopt the latter approach. It is convenient to rewrite ^ as: 

^"+1 =T3^ + TT^-^"' (2) 



yri+l = X. 



ni 



where a„ := \{zn + a), /3„ := ^(— z„ + a). Then we can regard ([2]) as a mapping 
defined on the domain F^ x F^. To resolve the indeterminacy at Xn = ±1, 
we apply the theory of the state of initial conditions developed by Sakai [5]. 
First we extend the domain to PF^ x PF^, and then blow it up at four points 
(x, y) = (±1, cxd), (cxd, ±1) to obtain the space of initial conditions: 

" ■— -^(1,00) '-'-^(-1,00) '-J -^(oo,!) ^-^(00,-1) ' W 

where -^^"00": ^^ ^^^ space obtained from the two dimensional affine space A^ by 
blowing up twice as 

-^("L):={((^-l>2/~').Ki:^i],["i:'^i]) 

?7i(a; - 1) = S,iy-\ (^i + a„??i)z;i = r/i(l - x)ui } C A^ x P x P, 



where [a : b] denotes a set of homogeneous coordinates for P-'^ . Similarly, 

772(2; + 1) == 6y"\ (-6 + Pnm)'"2 = '72(1 + X)U2 |, 

-^S,i)^={((^"''y-l)'K3:%],N:f3]) 

6(y - 1) = %2^~\ (% + a«C3)t'3 = 6(1 - y)u3. |, 



^ 



(") 



-1) •== {((^ ^2^ + 1)'K4 : '74], K : ^4]) 

^(y + 1) == rjiX^^, (-774 + ^„6)7^3 = 6(1 + y)u4 |- 



The bi-rational map ^ is extended to the bijection 0„ : QS"^^ -^ J7("+i) which 
decomposes as 0„ :— i„ o a)„. Here (,„ is a natural isomorphism which gives 
0(") = fi("+^) , that is, on Avl^ for instance, t„ is expressed as 

^ {{x ~ 1, jy-i), [e - ^2 . 77 : r,], [7. : 7;]) £ ^[^^^r 
The automorphism Con on fi'"-' is induced from ^ and gives the mapping 



A 



(n) 



l(") 



l(") 



,(n) 



l(") 



,(n) 



(1,00) ^ -^(004)' -^(cx),!) -^(-1,00)' •^(-1,00) ^ -^(00,-1)' •^(00,-1) ' •^(1,00)- 



^yl 



(«) 



Under the map -4["_^^) -J> -4(^_i), 



x = 1 



K 



(00,1) 



"3 



y 



/3. 



■"3, 



£;- 



(1,00) 



-^ E- 



(00,1) 



K 



(1,00) 



y' = i 



Ki : -^71] = [a«C3 + ^73 : 6] 



Vl 



2 ' 



where (x, y) £ -^(i 00) ' (^'' y') ^ "^(oo 1) ' -^i" ^^"^ -^2 ^'"^ ^he exceptional curves in 



i(") 



^p obtained by the first blowing up and the second blowing up respectively at 

\oo4) ^ -^(-1,00)' 



the point p S {(±1, cxd), (00, ±1)}- Similarly under the map A/^ -^^ -^ A/j^ 



El 



(00,1) 



El 



(-1,00) 



[6 : %] = [112 ■■ iPn - an)m - 6], 



K 



(00,1) 



^. 



(-1,00) 



[M3 : W3J 



-/3„7i2 : a„W2] 



The mapping on the other points are defined in a similar manner. Note that 
ujn is well-defined in the case a„ = or /?„ = 0. In fact, for «„ = 0, -E. 



(1,00) 



p(°oa) 



and E2 '' can be identified with the lines x = \ and y — 1 respectively. 



n" 



fi"" 



(0,0) J (»,») j 



(0,oo) (m,0) 



(1,0) 



(0,2) 



(1,2) 



(1,1) 



E 



(«,1),[1:2],[1:1] -►(2,»),[1:2],[1:2]- 
(oo,l),[l:2],[l:2]-^(2,»),[l:2],[l:l]- 
(«,1),[1:2],[0:1]->(2,»),[1:2],[0:1]- 

(l,co),[l:2],[0:l]->-(«,2),[l:2],[0:l]< 
(l,co),[l:2],[l:2]<-(«,2),[l:2],[l:l] < 
(1,»),[1:2],[I:1]<-(«,2),[1:2],[1:2]< 

(2,1) 

(0,1) ► (2,0) 



'1 
(2,2) 



(l,co),[l:0] 
(",2),[0:1] 



(»,1),[0:1] 
-(2,oo).[l:0] 



(»,1),[1:2],[1:0] 


(2,»),[1:2],[1:0] 
i 


(1,»),[0:I] 


(",2),[1:0] 

i 


(",2),[1:1] 


(1,»),[1:1] 
i 


(2,»),[1:1] 


(»,1),[1:1] 


(■»,1),[1:0] 


i 
(2,"),[0:1] 




i 


(1,»),[1:2],[1:0] 


(»,2),[1:2],[1:0] 



Figure 1: The orbit decomposition of the space of initial conditions fl^"'> and 
the reduced one 17^"' for p = r = 3. 



Therefore we have found that, through the construction of the space of initial 
conditions, the dPn equation can be well-defined over finite fields. However 
there are some unnecessary elements in the space of initial conditions when we 
consider a finite field, because we are working on a discrete topology and do 
not need continuity of the map. Let f2*^"^ be the space of initial conditions 
and |fi'^"^| be the number of elements of it. For the dPn equation, we obtain 
|f2(") I = (r + 1)2 - 4 + 4(r + 1) - 4 + 4(r + 1) ^ r^ + lOr + 1, since PF^ contains 
r + 1 elements. However an exceptional curve E^ is transferred to another 
exceptional curve Ef , and [1 : 0] S E2 to [1 : 0] £ i?| or to a point in Ef . 
Hence we can reduce the space of initial conditions f]'"^ to the minimal space of 
initial conditions $7^^"^ which is the minimal subset of fi'^"^ including PF^ x PF,., 
closed under the time evolution. By subtracting unnecessary elements we find 



|f7(")| = (r + l)2 -4 + 4(r + 1) 
following proposition: 



6r — 3. In summary, we obtain the 



Proposition 2.1 The domain of the dPu equation overVr can he extended to 
the minimal domain il^"^ on which the time evolution at time step n is well 



defined. Moreover |r2'"' 



6r — 3. 



In figure [1] we show a schematic diagram of the map tj„ on il^"\ and its 
restriction map a;„ :— CJn\n(") on fi*-"' with r = 3, ao = 1 and /?o — 2. We can 
also say that the figure[T]is a diagram for the autonomous version of the equation 



© when (5 = 0. In the case of r = 3, we have |f](")| = 40 and |fi(")| = 24. 

The above approach is equally valid for the other discrete Painleve equations 
and we can define them over finite fields by constructing isomorphisms on the 
spaces of initial conditions. Thus we conclude that a discrete Painleve equation 
can be well defined over a finite field by redefining the initial domain properly. 
In the next section, we show another extension of the space of initial conditions: 
we extend it to Z,, x Zjn. 



^p 



3 The discrete dynamical systems over a local 
field and its reduction modulo a prime (re- 
view) 

3.1 Almost good reduction 



, u 



Let p be a prime number and for each x ^ Q {x ^ 0) write x = p"p(^>— where 

V . 
Vp{x), u, w e Z and u and v are coprime integers neither of which is divisible by 

p. The p-adic norm \x\p is defined as \x\p = p'^p'^^y (|0|p — 0.) The local field 

Qp is a completion of Q with respect to the p-adic norm. It is called the field 

of p-adic numbers and its subring Zp := {x e Qp| |a;|p < 1} is called the ring of 

p-adic integers [17j . The p-adic norm satisfies a non- Archimedean (ultrametric) 

triangle inequality |a;-|-y|p < max[|a::|p, \y\p\. Let p = pZp = {x £ 'Lp\ Vp{x) > 1} 

be the maximal ideal of Zp. We define the reduction of x modulo p as 

X '.Ijp 3 X i-^ X :~ X mod p G Zp/p = Fp. 

Note that the reduction is a ring homomorphism. The reduction is generalised 

toQ^: 



i;3. = pS.(«ez;)^{^ [fcJo) 



which is no longer homomorphic. Note that for a; G Qp , A: > is equivalent to 
a; G Zp . Let us consider a rational system (p consisting of two rational functions 
of X and y: 



4'ix,y) = ((/)i(a;,2/),02(a;,j/)) G (Zp(x,j/)) : Qp 
its reduction 



' '% 



0(a;,y) = {(i>iix,y),(f)2ix,y)) G (Fp(x,?/))2 

is defined as the system whose coefficients are all reduced. The rational system is 
said to have a good reduction (modulo p on the domain V) if we have (t>{x, y) = 
(j){x,y) for any {x^y) G T) ^. We have defined a generalised notion in our 
previous letter and have explained its usefulness; 



^With the method in the next section, we can only deal with systems over the field F i 
To apply it to those over the field Fpm (m > 1), we have to use the field extension of Qp. 



Definition 3.1 ([!]) A (non- autonomous) rational system (j)n: Qp -^ Qp {n G 
Z) has an almost good reduction modulo p on the domain D C Z^, if there exists 
a positive integer rUp-n for any p = {x,y) G V and time step n such that 

<l>n'"" {x, y) = 0™"^" (x, y), (4) 

where 0™ := 0„+„j_i o 0„+„j_2 o • • • o 0„. 

Let us first review some of the findings in [1] in order to see the significance of 
the notion of almost good reduction. Let us consider the mapping ^-yi 

axn + 1 

where a G {1, 2, • • • ,p — 1} and 7 G Z>o are parameters. The map ([5]) is known 
to be integrable if and only if 7 = 0, 1, 2. Note that when 7 = 0, 1, 2, ^ is 
an autonomous version of the g-discrete Painleve I equation and therefore is 
integrable. We also note that when 7 = 2, ([5]) belongs to the QRT family and 
is integrable in the sense that it has a conserved quantity. 

Let V be the domain {{x, y) £ "ii \x ^ Q,y j^ 0}, then clearly 

*2(a;„,2/„) = *2(in,yn) for Xn ^0, jjn ^ 0. 

For (x„,y„) G V with a;„ = and j/„ ^ 0, we find that ^2(^1' = Oii/n) is not 
defined for k = 1,2, however it is defined if fc = 3 and we have 

mxn,yn) ^ min -- 0,y„) - (^^,0 

Finally for Xn ^ yn ~ 0, we find that ^2(2^": Vn) is not defined for k = 1,2, .., 7, 
however ,.-^_^ 

*|(a;„,y„) = *i(i„ =0,y„ = 0) = (0,0). 

Hence the map ^2 has almost good reduction modulo p on V. Note that, in the 
case 7 = 2 and a = 0, if we take 

/2fc := X2kX2k-l, flk-l := {x2k-lX2k-2) 

([5]) turns into the trivial linear mapping fn+i — fn which has apparently good 
reduction modulo p. In a similar manner, we find that ^^ (7 = 0, 1) also has 
almost good reduction modulo p on T). On the other hand, for 7 > 3 and 
in = 0, we easily find that 

^A:GZ>o, *;^(a;„,2/„)7^*;^(i„ = 0,y„), 

since the order of p diverges as we iterate the mapping. Thus we have proved 
the following proposition: 

Proposition 3.2 The rational mapping ([Sj has an almost good reduction mod- 
ulo p only for 7 = 0, 1, 2. 

Note that having an almost good reduction is equivalent to the integrability of 
the equation in these examples. 
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3.2 The dPn equation modulo a prime and its special so- 
lutions 

Now let us examine the dPn ([2]) over Qp. We suppose that p > 3, and redefine 
the coefficients a„ and /3„ so that they are periodic with period p: 

{iS + Zf) + a + Uap) {-iS - zo + a + npp) 

(meZ, ie {0,1,2,--- ,p-l}), 

where the integer n^ (np) is chosen such that £ {a^lfJo (0 G {/^ilfJo )■ As a 

result, we have a„ = "'^+f+° , /3„ = ^"'^-^°+° and |a„|p, |;3„|p £ {0, 1} for any 
integer n. 

Proposition 3.3 Under the above assumptions, the dPu equation has an al- 
most good reduction modulo p on D := {{x, y) £ Z^ \x ^ il}- 

proof We put {Xn+l,yn+l) = 4'n{Xn,yn) = i4>n{Xn,yn),(f>n{Xn,y7i)]- 

When Xn y^ ±1, we have 

dn , /3n . _ 

a^n+l = _ . + ——: Vn, yn+1 = Xn- 

Hence 0„(a;„,j/„) = 4>n{xn,yn)- 

When Xn — 1, we can write x„ = 1 +^^^6 (fc £ Z+, |e|p = 1). We have to 
consider four caseo 
(i) For a„ = 0, 



in+l = 4>n^ {Xn.Vn) = ( ^ ) ~ J/"' 



Hence we have (/)„(x„,y„) = 4>n(xn,Vn)- 
(ii) In the case a„ 7^ and /3n+2 7^ 0, 



_ (Qn -/3n)(l+ep*^) +Q _ _ 2q„ + {an - Pn)ep^ 
ep'^(2-\-ep'^) ep''[2 + ep'^) 

a^ + polynomial of 0{p) 
a^ + polynomial of 0{p) ' 
_ {2a„y„ + 2Sl3n+i + (2 - (5)a}a^ + polynomial of 0{p) 
2/3„+2a^ + polynomial of 0{p) 

Thus we have 

_ 23„y„ + 2(5/3„+i + {2 - d)a , _ 

^n+3 — —3 1 2/n+3 — ~-L: 

2pn+2 



■^Precisely speaking, there are some special cases for p = 3, 5 where we have to consider the 
fact a„+p = On or l3„+p = fin during the iteration process. We can prove by straightforward 
calculations that similar results hold for these exceptional cases. 



and (/)3(x„,y„) = (j)l{xn,yn)- 

(iii) In the case «« ^ 0, /3n+2 = and a ^ —S, we have to calculate up to Xn+5- 

After a lengthy calculation we find 



Xn+4 = f/*^ (l,2/«) = 1, and J„+5 = 05 (l,y„) = -- 

and we obtain 0^(x„, y„) = (/>^(i„, y„). 

(iv) Finally, in the case a„ 7^ 0, /3„+2 = and a — —S we have to calculate up 

to x„+7. The result is 

a;n+6 = 0„ (1, yn) = -1, and Xn+7 = K (1: yn) = ^ , 

and we obtain (/)^ (x„ , j/„) = (^^ (£„ , j/n) ■ Hence we have proved that the dPn equation 
has almost good reduction modulo p at a;„ = 1. We can proceed in the case 
Xn = — 1 in an exactly similar manner. □ 

From this proposition, the evolution of the dPn equation ([T]) over PFp can 
be constructed from the initial values u„_i and u„. 

Now we consider special solutions to (IT|) over PFp. For the dPu equation 
over C, rational function solutions have already been obtained [TO]. Let A'' 



2(JV+1) 
A' 



be a positive integer and A ^ be a constant. Suppose that a — — j- 
6 = zo = ^ 




(fee z<o). 



and 



T^ := det 



^N+l-2i+jW 



i<iJ<N 



(6) 



Then a rational function solution of the dPn equation is given by 



u„ = 



n+l _n-l 

' n+i'n 



1. (7) 



'JV+l'JV 



If we deal with the terms in (|6]) and ([7]) by arithmetic operations over Fp, we 
encounter terms such as - or - and ^ is not well-defined. However, from 
proposition 13. 3[ we find that ^ gives a solution to the dPn equation over PF^ 
by the reduction from Q (c Qp), as long as the solution avoids the points 
(a„ = 0, Un — 1) and (/3„ =0, Un — —1), which is equivalent to the solution 
satisfying 

N+l N-l 
^-N-l^-N-3 , n iV + 1 N , r, ic,\ 

^N+\ ^N ^ L)' N W- ^ ^' \^) 

~N+l~N 

where the superscripts are considered modulo p. Note that r]^ = '''w for all 
integers TV and n. In the table below, we give several rational solutions to the 



dPji equation with N = 3 and A = 1 over PF^ for q = 3,5,7 and 11. We see 
that the period of the solution is p. 



p 


-JV-l -JV-3 


JV + l N-l 
N + l JV 

N + 1 N 




111, M2, «3,"4. "5, «6,M7. «S, "9,"10: ■ ■ • 


3 


OO 


OO 


1, 2, OO , 1, 2, OO, 1, 2, OO, 1, 2, OO, 1, 2, OO, . . . 
period 3 


5 


OO 


4 


4, 2, 3, 1, OO, 4, 2, 3, 1, oo, 4, 2, 3, 1, oo, 4, . . . 
period 5 


7 


OO 





1, oo, 6, 5, 1, oo, 6, 1, oo, 6, 5, 1, oo, 6, 1, oo, . . . 


period 7 


11 





7 


oo, 1, 6, 1, oo, 10, oo, 1, 0, 2, 10, oo, 1, 6, 1, . . . 


period 11 



We see from the case of p = 11 that we may have an appropriate solution 
even if the condition ([8]) is not satisfied, although this is not always true. 
The dPii equation has linearized solutions also for 6 — 2a [20]. With our new 
method, we can obtain the corresponding solutions without difficulty. 



4 The j9-adic analogue of the singularity confine- 
ment test 

The above approach is closely related to the singularity confinement method 
which is an effective test to judge the integrability of the given equations [T5] . 
In the proof of the Proposition l3.3[ we have taken a;„ = 1 + ep'^ and have shown 
that the limit 

Um yXn-\-m7 '^n-^m-{'l } 
|ep'=|p^O 

is well defined for some positive integer m. Here ep^ (k > 0, \e\p = 1) is an 
alternative in Qp for the infinitesimal parameter e in the singularity confinement 
test in C Note that p*^ {k > 0) is a 'small' number in terms of the p-adic 
metric. From this observation and propositions 13. 2[ 13. 3[ we postulate that 
having almost good reduction in arithmetic mappings is similar to passing the 
singularity confinement test. 



5 The g-discrete Painleve equations over a finite 
field 

In this section we study the g-discrete analogues of the Painleve equations 
{qPi and gPn equations). 
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5.1 The qPi equation 

One of the forms of the qPi equations is as fohows: 

_ ag"x„ + b 

•^n+l-^n — 1 — 2 ' V / 

where a and b are parameters. We rewrite ^ for our convenience as 

aq"xn + b 

'^n:{ "^"^^ 2 xlVr. ' (10) 

yn+1 — '^n- 

Similarly to the dPn equation, we can prove the following proposition: 

Proposition 5.1 Suppose that a,b,q are integers not divisible by p, then the 
mapping (jlOp has an almost good reduction modulo p on the domain P := 
{ix,y)eZl \x^O,y^O}. 

proof Let (a:„+i,y„+i) = ^n{xn,yn)- Just like we have done before, we 
have only to examine the cases i„ = and j/n = 0. We use the abbreviation 
q = q,d = a,b = b for simplicity. By direct computation we obtain; 
(i) If Xn — and y„ ^ 0, then 

~ / 52 

(ii) If y„ = and x„ ^ 0, then 

$5(x„,2/„) = $5(i„,0)= fo,||- 
(iii) If in = and y„ = 0, then 

<i>«'(^J/n) = ^(0,0) = (0,0). 

5.2 The Q'Pii equation 

The qPii equation is the following g-discrete equation: 

iz{qr)z{r) + l)(z(r)z(g- V) + 1) = fl^lli^. (n) 

T — z{t) 

where a and g are parameters [21j . It is also convenient to rewrite (jlip as 

_ Q(g"To)^Xn - (g"ro - Xn)(l + XnZ/n) 

*n : ■( "^^ a;„(g"To -a;„)(x„y„ + 1) ' (12) 

where r = g'Vo. Just like the qPi equation, we can prove that qPn has an 
almost good reduction: 
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Proposition 5.2 Suppose that a,q,TQ are integers not divisible by p, then the 
mapping (|12p has an almost good reduction modulo p on the domain V := 
{{x,y)GZl \x^O,x^q^To (nGZ),xy + 1^0}. 

proof Let {xn+i,yn+i) = '^nixniVn)- Just like we have done before, we have 
only to examine the cases £„ = 0, g"To and —Vn^- We obtain; 
(i) If in = and -1 + g^ — aq'^r^ + q^r"^ — q^ryn t^ 0, 



^l{xn,yn) = $^(i« = 0,y„) 



1 — g + aq T — q'^r — aq t^ + q Ty„ g 

q T 



(7^t(— 1 + q'^ — aq^T^ + q^r"^ — q'^TJjn) 
(ii) If in = and — 1 + g^ - ag^r^ + q^r^ - q^ryn = 0, 

' 1 — g^ + q^T^ — aq^r"^ 



*n(a;«,yn) = *^(i„ = 0,y„) 
(iii) If in = r and 1 + ry^ 7^ 0, 



q'^T 



-,o 



^n{Xn,yn) = *^(in =T,yn) 

^ n-q^ + {a + q- aq^)q^T'^ + (1 - g^)r^ + (1 - aq)q^T^y 
(iv) If in = T and 1 + Tyn = 0, 

(v) If i„y„ + 1 = 0, 



,0 



19 "^ ' 



-aq'V^ 



-^ 12 4~ 

Thus we complete the proof. □ 

From this proposition we can define the time evolution of the qVn equation 
explicitly just like the dPn equation in the previous section. 

Next we consider special solutions for qPn equation ([TT|) over PFp. In [22] 
it has been proved that (ITT]) over C with a = g^jv+i ^^ g ^-j jg solved by the 
functions given by 



0(^)(r) 



g^+igW(gr)5(^+i)(r) 



w(<7 



-1+2J-1, 



l<ij<Af 



5^-^(r) = 1 



w{q 



i-2i, 



l<ij<-Ar 



(^>0), 
(^ < 0), 

(7V>0), 
(iV = 0), 
(A^<0), 



(13) 



(14) 
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where wIt) is a solution of the g-discrete Airy equation: 
w{qT) — tw{t) + w{q^ r) — 0. 



(15) 



As in the case of the dPn equation, we can obtain the corresponding solutions 
to P^ over PFp by reduction modulo p according to the proposition 15.21 For 
that purpose, we have only to solve (J15p over Qp. By elementary computation 
we obtain: 

wiq^'^^To) = ciP„(ro;q) + coF„-i(qTo; g), (16) 

where cq, ci are arbitrary constants and Pn{x;q) is defined by the tridiagonal 
determinant: 



Pn{x;q) 



qx —1 

— 1 q'^x 







q''-^x 





-1 



The function Pn{x; q) is the polynomial of nth order in x, 



V2] 



P„(x;(7)= V(-l)V„;fc(q)x" 



-2fc 



fc=0 



where an;k{_<l) are polynomials in g. If we let i <^ j denotes i < j — 1, and 
c{ji,J2, ■■■,3k) ■= Z]r=i(2jV + 1), then, we have 

a„.fc = Y, qn(n+l)/2-c(n,n,-,0u)_ 

l<il«i2«<Cjfc<n-l 

Therefore the solution of qPn equation over PFp is obtained by reduction modulo 
p from (HH), (dH) and ^ over Q or Qp. 

We can also prove that a type of g-discrete Painleve III, IV and V equations 
have an almost good reduction on some appropriate domains. These facts jus- 
tifies that almost good reduction can be a criterion for integrability of systems 
over finite fields. Proofs are fairly straightforward and will be reported in the 
future article [23] . 

6 The discrete KdV equation 

In section[21 we have successfully determined the time evolution of the dPn equation 
through the construction of the space of initial conditions of the dPn equation 
by blowing-up twice at each of the singular points so that the mapping becomes 
bijective. However, for a general nonlinear equation, explicit construction of 
the space of initial conditions over a finite field is not so straightforward (for 
example see |16) ) and it will not help us to obtain the explicit solutions. 
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In this section we study the sohton equations evolving as a two-dimensional 
lattice over finite fields by following the discussions made in |3J. For example, 
let us consider the discrete KdV equation 

' ^ + T^K+^--J.+i)-0- (17) 






over a finite field F^ where r = p™, p is a prime number and m £ Z+. Here 
n, t £ Z and (5 is a parameter. If we put 

— -^ {I + S)^^ - 6xi 

we obtain equivalent coupled equations 

(1 + 5)yi 



Vn+l — 



1 + Sx^y*^ ' 
(1 + Sxl,yl)xi 



(18) 



1 



Clearly (fT8|) does not determine the time evolution when l + 6x*„yl^ = 0. Over a 
field of characteristic such as C, the time evolution of (x^, yj^) will not hit this 
exceptional line for generic initial conditions, but on the contrary, the evolution 
comes to this exceptional line in many cases over a finite field as a division by 
appears. The mapping, (a;*^, yj,) n- (x^+^, 2/^^2)1 i^ lifted to an automorphism of 
the surface X, where X is obtained from P^ x P^ by blowing up twice at (0, 00) 
and (00,0) respectively: 



X = Am no^ U A 



(0,00) U 'f^(cx),0)) 



C A^ X pi X 



A(o^oo) := |((a;,y ^) ,[^-'n],[u-v]) x-q = y ^^,r]u = y ^{i] + S^)vj 

A(oo,o) ■■={{{x-\y),[^:Tj],[w:z]) x'^rj^y^, {?] + d^)u = yj^v} C A^ x P^ > 

To define the time evolution of the system with N lattice points from ([TS)) . 
however, we have to consider the mapping 

{yi;xi,x2, ■■■,xpf) I > [Xi ,X2 ,...,x^ ;y7v+i)- 

Since there seems no reasonable decomposition of X into a direct product of two 
independent spaces, successive use of (fT8)) becomes impossible. Note that if we 
blow down X to P^ x P^ , the information of the initial values is lost in general. 
If we intend to construct an automorphism of a space of initial conditions, 
it will be inevitable to start from P^+i and blow-up to some huge manifold, 
which is beyond the scope of the present paper. This difficulty seems to be 
one of the reasons why the singularity confinement method has not been used 
for construction of integrable partial difference equations or judgement for their 
integrability, though some attempts have been proposed in the bilinear form 
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ii 



PB] . There should be so many exceptional hyperplanes in the space of initial 
conditions (if it does exist), and it is practically impossible to check all the 
"singular" patterns in the naive extension of the singularity confinement test. 
On the other hand, when we fix the initial condition for a partial difference 
equation, the number of singular patterns is restricted in general and we have 
only to enlarge the domain so that the mapping becomes well-defined. This is 
the strategy that we will adopt in this section. 

First we explain how the indeterminate values appear through the time evo- 
lution by examining the discrete KdV equation ([T^ over F7 :— {0, 1, 2, 3, 4, 5, 6}. 
If we take S ~ 1, HH]) turns into 



vU 



1 + <yi ' 

(1 + x\,yi)x\^ 



2 
Suppose that x\ — &, X2 = 5, yj = 2, y\ = 2, then we have 

a;} = ^^3, 2/2 = f ^4 mod 7. 

With further calculation we have 

2_4^4 i_21^ ,_8^1 

^1-7-Q, y2- 2 -'J' ^2-21-0- 



4 1 

Since - and - are not defined over F7, we now extend F7 to PF7 and take 

j 

— = 00 for i e {1, 2, 3, 4, 5, 6}. However, at the next time step, we have 

2_ 2-0 ^_(I + cx)-0)-cx) 

"^2= l-t-oo-O' ^3= ^ 

and reach a deadlock. 

Therefore we try the following two procedures: [I] we keep 5 as a, parameter 
for the same initial condition, and obtain as a system over ¥j{5), 

1 _ 2(1 + 5) o_ 6{l + h5) 
^1 " 1 + 5(5 ' ^2 " 1 + 5 ' 
1 _ 6(l-h(5)(l + 5(5) 1 _ 2(1 + 25 + 4(52) ^ _ 2{l + 5){l + b5) 



„2 



1 + 3(5 + 3^2 ' ^2 ^^ ^ 5^^2 ' ■ i 1 + 2(5 + 4(52 

4(l + (5)(2 + (5)(3 + 2(5) ^ _ 2(5 + 5,5 + 2J2) 



2 



"^2 (l + 5(5)(5 + 5(5 + 2(52)' ^^ (2 + (5) 

[II] Then we put (5 = 1 to have a system over PF7 as 



.} = 


= 3, 


2/2° = 


4, 


x\ 


= 


9 


120 




1 




24 


X2-- 


72 


= 4, 


vl 


— 


¥ 



72 _ 1 14_ 2 24_ 

= 00, Vn = = U, Xt — = CXD, 
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y° = 2 

xl = 3 x\ = ca ^3=0 

yi = 2 

X? = 4 



^^ y2° = 4 ^U y3° = ^U 



^^ yi = o ^^ y3^ = 5 ^^ 



Figure 2: An example of the time evolution of the coupled discrete KdV equation 
(fTS]) over PF7 where 5^1. 



Thus all the values are uniquely determined over PF7. Figure [5] show a time 
evolution pattern of the discrete KdV equation p^ over PF7 for the initial 
conditions x\ = 6, x^ = 5, xl =A, xl^ 3, x° = 2 (j > 5) and y\ = 2{t> 0). 
This example suggests that the equation ([18| should be understood as evolv- 
ing over the field Fj.((5), the rational function field with indeterminate 5 over F^. 
To obtain the time evolution pattern over PF^, we have to substitute 5 with a 
suitable value 5q € F^ (<5o = 1 in the example above). This substitution can be 
expressed as the following reduction map: 

( (s > 0), 

F.(^)x^PF„: (^-^or§^- Z ^'^'^' 

f{S-So) g{0) 

[m ^'^'^' 

where s G Z, f{h), g{h) G ¥r[h] are co-prime polynomials and /(O) ^ 0, 5(0) ^ 0. 
With this prescription, we know that 0/0 does not appear and we can uniquely 
determine the time evolution for generic initial conditions defined over F^. 

6.1 Soliton solutions of the discrete KdV equation over 
finite fields 

We consider the A^-soliton solutions to ([T7| over F^. As mentioned in the intro- 
duction, the A^-soliton solution is given as 



l-hY f h+S 



X* 




-* +1-^' ' 


< 


:= det 5ij 





h + lj-l 




16 







1 



9 10 oc 



■BTT 



Figure 3: The one-soliton solution of the discrete KdV equation ([T7|) over R 
(left) and PFn (right) where ,5 = 7, 71 = 2, /i = 9. Elements of PFn are 
represented on the following grayscale: from (white) to 10 (gray) and 00 
(black). 



where 7^, h [i = 1,2,...,N) are arbitrary parameters but li ^ Ij for i j^ j. 
When li, 7i are chosen in F^, x*„ becomes a rational function in Fr((5). Hence 
we obtain soliton solutions over PF^ by substituting S with a value in F^. The 
figures [3] and U show one and two soliton solutions for the discrete KdV equation 
P7)) over the finite fields PFn and PFig. The corresponding time evolutionary 
patterns on the field M are also presented for comparison. 

The soliton solution of the discrete KdV equation over the finite field F^ has 
a period r — 1 since we have a^~-^ = 1 for all a £ F^ . 



6.2 The discrete KdV equation over Qp 

We can also define the time evolution of the discrete KdV equation over the 
field Fp by considering the system over the field of p-adic numbers Qp instead of 
Fp((5), in the same way as we have done to the discrete Painleve equations. We 
just have to compute the evolution of (fT8|) or its soliton solutions over Qp and 
then reduce the results to Fp. (To deal with the equation over Fpm (to > 2), we 
require the field extension of Qp.) 



7 Concluding remarks 



In this article we investigated the discrete Painleve equations and discrete KdV 
equation over finite fields. To avoid indeterminacy, we examined two approaches. 
One is to extend the domain by blowing up at indeterminate points. According 
to the theory of the space of initial conditions, this approach is possible for 
all the discrete Painleve equations. An interesting point is that the space of 
initial conditions over a finite field can be reduced to a minimal domain because 
of the discrete topology of the finite field. The other is to define the system 
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ruhx 



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 



Figure 4: The two-solitoii solution of the discrete KdV equation (|17p over R 
(left) and PFig (right) where S = 8, 71 = 15, h = 2, 72 = 9, h = 4. Elements 
of PF19 are represented on the following grayscale: from (white) to 18 (gray) 
and cx) (black). It is difficult to see the interaction of solitons over PFig. 



over the larger field (Qj, in the case of dPn , qPi and gPn , or Wp{6) in the 
case of discrete KdV equation), and then reduce the results to finite fields. 
In particular the systems over Qp attracted our attention. We defined a new 
integrability test {almost good reduction), which is an arithmetic analogue of 
the singularity confinement test, and proved that several types of the discrete 
Painleve equations have this property. Thanks to this property, not only the 
time evolution of the discrete Painleve equations can be well defined, but also a 
solution over Q (or Qp) can be directly transferred to a solution over PFp. This 
approach is equally valid in other discrete Painleve equations and we expect 
that the same is true for its generalisations [M]. Furthermore, this 'almost good 
reduction' criterion is expected to be applied to finding higher order integrable 
mappings in arithmetic dynamics. In the last section, we studied the discrete 
KdV equation and showed that a similar approach is also useful in defining 
the discrete partial difference equations such as soliton equations over finite 
fields and in obtaining the soliton solutions of these equations. One of the 
future problems is to investigate the property of solitary waves of discrete KdV 
equation and other soliton equations such as discrete KP equation and discrete 
nonlinear Schrodinger equation over the finite field. Solving the initial value 
problems for these equations over the field of p-adic numbers, the field of rational 
functions and the finite fields is also an important problem to be studied. 
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